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The mounting of either a tractor or pusher propeller near 
a streamlined wing not only affects the efficiency of the pro- 
peller, hut also produces variations in the flow which, being 
always unsymmetrical, may subject the propeller to additional 
aerodynamic stresses. The purpose of the present investigation 
was therefore to determine the variations in the velocity of the 
air when it was obliged to flow about a wing profile. Of espe- 
cial interest was the determination of the velocity as compared 
with that of the undisturbed flow and of the velocity drop in 
successive planes perpendicular to the flow direction at increas- 
ing distances from the wing profile. . On the basis of these re- 
sults, it is then possible to obtain an approximate picture of 
the asymmetry of the loading of the propeller disk at different 
positions of the propeller with respect to the wing profile. 

In order to enable a theoretical calculation of the flow 

relations about a profile, we must choose profiles susceptible 

of hydrodynamic and mathematical investigation, like the Joukows- 
*"Untersuchungen uber die Veranderung der Windgeschwindigkeit bei 
Umstromung eines .Flugelprof iles. " From Zeitschrift fur Flugtech- 
nik und Motorluf tschif f ahrt, July 15, 1939, pp. 341-347. 
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ky profiles. The possible variants in the construction of such 
a profile were so correlated that the resulting chord, thickness 
and camber approximately agreed with the corresponding dimen- 
sions of the profiles in actual use* 

By taking such a profile as the basis, it is possible, with 
the aid of mathematical formulas employed in conformal transfor- 
mation for the development of wing profiles from a circle, to 
transfer the relatively' easily determined flow relations with 
respect to the course and the velocity to the corresponding 
points of the field of flow. about the J-profile, thus rendering 
it possible to determine the theoretical flow relations in the 
vicinity of the profile to be investigated. It is not possible 
by theoretical investigation to determine the effect of the 
boundary zone in the immediate vicinity of the profile, where 
laminar flow conditions prevail which are not directly suscepti- 
ble of theoretical calculation. The boundary zone in question 
is so small, however, that it can be disregarded in practice. 

The selected J-profile is the Gottingen profile No. 433. 
It has the characteristic parameter for the camber, nsunely, the 
approximate ratio of the rise of a circular arc, plotted as a 
curved profile axis, to the corresponding half-chord f/l =0.1. 
For the thickness, the approximate ratio of the maximum profile 
thickness to the profile chord d/l = 0.15. This profile was 
constructed from these values, as shown in Figure 1. 
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The angle of attack then had to be determined for this pro- 
file, since the effect of the flow velocity depends on this an- 
gle. It appeared to be the most important to investigate the 
flow relations at sea level, for which an angle of attack of 0°, 
or somewhat less, has been found to- be the best. The adoption 
of this angle of attack gives the same angle of attack for the 
generating fundamental circle and establishes the position of 
the points of. full dynamic pressure on which the calculation of 
the circulation and the course of the. {^symmetrical flow. about 
the circular profile is based* The details of this method are 
fully explained in the later writings of E. Trefftz and Oskar 
Schrenk. Hence only a. brief summary will be given here* The 
un symmetrical flow about the circular profile, which corre- 
sponds to the determining flow for the wing profile, may be con- 
sidered as the superposing of a symmetrical flow about the circu- 
lar profile and a circulatory flow about the same. The general 
potential equation of the two-dimensional flow about a circular 
profile is written as follows: 

For the symmetrical flow, 

f(z) = - v.' (z + 5*1) (i) 

For the circulatory flow, 

f(z) = - t jpft m z (3) 

Z represents the complex distance of any point in the Z-plane 
from the origin of the coordinates and is therefore 
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The introduction of complex coordinates is necessitated by 

the principle of conformal or orthomorphic transformation. 
Vp, represents the flow velocity outside the disturbed region 

around the cylinder. 
R x represents the radius of. the fundamental or basic circle 

(R x = X x 2 + Y x 2 ). 
T represents the circulation corresponding to the location. of 

the center of dynamic pressure on the fundamental circle. 

With the introduction of the coordinates X + iY in the 
above functional equations and the division into real and imagi- 
nary parts, the former represents the potential and the latter 
the stream function. We obtain for them, with a few transfor- 
mations in the second case, the following equations: 

*(X,Y) -▼•(!- £\ ]i ) (3) 

• + (X,Y) = - m/ X 3 + Y 2 (4) 

By putting each of the two functional equations equal to a 
constant C x or C 2 , we can determine the course of a streamline, 
if we solve the equations for X and Y and calculate, for the 
assumed value of one of the unknowns, the corresponding value of 
the other unknown. Variations in the value of C thus enable 
the determination of a whole group of streamlines. Equation (4) 
requires only a knowledge of the magnitude of T, or V can be 



N.A.C.A. Technical Memorandum Nq. 575 5 

replaced by an expression containing only quantities which can 
be measured from the drawing. This result is attained by the 
following line of reasoning. 

By differentiating the potential functions, we obtain the 
velocity functions, which give for equations (l) ajLd (2), after 
suitable transformations, the following expressions for the ab- 
solute velocity: 



W 



2 V sin A 

00 



217 y x 3 + y 3 



for the parallel flow (5) 
for the circulatory flow (6) 



The angle A is here formed by a line parallel to the direc- 
tion of flow through the center of the fundamental circle and a 
line from the same point to the corresponding point of the field 
of flow with the coordinates X and Y. Equation (5)' is a limit- 
ing form and is fulfilled for all points of the fundamental cir- 
cle. It also holds good for the centers of dynamic pressure on 
the same circle. At this point, however, the velocity of the 
parallel flow must equal the velocity of the circulatory flow, 
since the velocity is zero at the points of full dynamic pressure. 
For these points we therefore have 

v r 

2 TL sin A = === _ = — ' (?) 

For this point the angle A has a definite value a th which 
can he measured on the drawing. Hence we have 
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T = 4 Tt R 1 V tt sin a th (8) 



or 



i(X,Y) = 2 R x V. sin In J I s + Y 2 (9) 



Both groups of curves for the circulatory flow and for the 
parallel flow can" now he calculated and plotted by transforming 
the equations (3) and (4)to the^form 



x = / - °* Y - Y R * (10) 




9 C * ( 11 ) 



x _ / e 2 Ri sm cc th _ y3 



The last equation has the form of the general equation for 
a circle. Hence, on putting the term Y 2 equal to zero, X 
gives the radius: Q 

X = e 2 Rl 83111 °* h • (12) 

The curves in Figure 2 are obtained from equations (10) and 
(18). Their calculated values are given in Tables I and II. 
Attention is called to the fact that, for the numerical calcula- 
tion, R x is taken as the unit of mass and accordingly equals 
1, and the values of I and Y are therefore multiples of the 
cylinder radius. 

The two sets of curves plotted from these numerical values 
correspond therefore to the parallel flow and to the circulatory 
flow, but not to the unsymmetrical flow, which must satisfy the 
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functional equation 

*(x,t) = v,, (y - - 2 R i sin a th m/?^ y 2 ) (is) 

The numerical solution of this equation is extremely diffi- 
cult, so that it is possible to calculate only individual points 
of irregularly staggered flow constants by the introduction of 
any convenient X and Y values. The streamlines can be easily 
determined, however, by the graphic combination of both members 
of the equation. Hence. the above-described method was adopted 
at the outset. Figure 2 also shows the new set of curves for the 
unsymmetrical flow, as derived from the two original sets. 



TABLE I 



/ Y 3 - C Y 2 - Y R. _ 
X = J ^—^ i; C < Y, R x = 1 



0, = o 


Y 

X 


0 1 

1 0 


Circumference of Circle 


C\ = 0.1 


Y 

X 


1.05 

0 


0.75 
0.768 


0.5 
1.0 


0.25 
1.305 


0.15 
1.725 


0.14 
1.86 


. 0.1 

CD 


0 X = 0.25 


Y 
X 


1.133 
0 


1.0 
0.576 


0.75 
0.969 


0.5 
1.32 


0.3 
2.43 


0.25 

00 




0 X = 0.4 


Y 
X 


1.22 

0 


1.0 

0*815 


0.75 
1.26 


0.5 

2.18 


0.4 

00 






Oj = 0.6 


Y 
X 


1.34 

0 


1.0 

1.225 


0.75 
2.04 


0.65 
3.54 


0.6 

00 






Ci = 0.8 


Y 
X 


1.478 

0' 


1.25 
1.105 


1.0 
2.0 


0.85 
4.04 


0.8 

00 






Ci = 1.0 


Y 
X 


1.62 

0 


1.50 
0.872 


1.25 
1.85 


| 1.0 . 

00 








Oj = '1.3 


Y 

X 


1.77 

0 


1.5 
1.66 


1.2 

00 










Ox = 1.4 


Y 
X 


1.92 

0 


1.75 
1.41 


1.6 

2.34 


1.4 

CO 








Ox = 1.6 


X 


2.08 

0 


2.0 
1.0 


1.6 

00 










Ox = 0.05 


Y 
X 


0.125 
1.284 
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TABLE II 



X = e 2 Rl Sin ^ sin c^ h = 0.1583} R x = 1 

log X = CL 0 l0g | ; log X = 0 2 1*3*1 
& ^ 2 sin 0L. h . ' . 3 . 





0 


0.05 


0.1 


0.15; 


0.2 


0.25 


0*3 


0*4 


0*5 


0*55 


Y 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


X 


1 


1.171 


1.371 


1.605 


1.881 


2.201 


2.580 


3.536 


4.849 


5*6.75 



Moreover, the course of the flow around the basic circular 
profile is well known. According to the principle of conformal 
transformation there corresponds to every point in the Z plane 
outside the basic circle a definite point within the auxiliary 
circle of the profile construction, whose location is determined 
by the following conditions. 

1. To. the line from the point 0 of the profile construc- 
tion: there corresponds a line mirrored with respect to the hori- 
zontal and passing through the zero point, which serves as a 
common axis of reference. 

2. The length of this line is f/Z, where Z denotes 
the length of the upper line. 

The bisection of the line between these two points gives 
the location! of the point of the z-plane corresponding to the 
outer point of the g-plane for which the profile is also con- 
structed. Strictly speaking, we ..obtain by this profile construc- 
tion all the measurements in the z-plane only half-scale, but 
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this is of no consequence in the present velocity calculation, 
since it is corrected by using the coefficient 2. The syste- 
matic transposition of a series of points of a streamline enables 
us to find the corresponding flow about the wing profile, in 
which connection attention may be called to the fact that every 
point found is correlated with the corresponding point on the 
circle and therefore has a flow velocity comparable with that of 
the corresponding point of origin* If the velocity of the un- 
symmetrical cylindrical flow is thus known, the velocity for the 
corresponding profile streamline can be easily found. The above 
method was used for a number of streamlines, the results being 
given in Figures 3 and 4. 

The equation for the flow velocity at any given point was 
previously developed. For the determination of the complex po- 
tential equation, we will start once more from the complex poten- 
tial equation. 

R 2 

f(Z) = - v. (z + -~) 

If we put Z = R e iA , that is, if we introduce polar coordinates 
with respect to the center of the circle, equation (14) is 
changed to 

R s \ 

» p - - t- (i-pinr) 

or 

lf p = . V a (1 - i cos 2 A + i m sin 3 A) (15) 
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if we put 




(16) 



If we put 



(18) 



c = 1 - m 



the absolute value • of • TIL. • is ■ calculated' with 



W_ =V / (l-- m cos'2 A) • +-m 'Sin 3 2 A • 





(17) 



In equation (16) 
R x denotes the radius of the basic circle; 

R denotes the length of the line from the center of the basic 
circle to the streamline point for which the velocity is 
to be calculated; 

A denotes the angle between the line R and the parallel to 
the flow direction through the center of the basic circle. 

If we measure the vertical distance h of the point on 
the streamline from, the parallel to the direction of flow, we 
obtain 



Hence all the requisite numerical values can be determined 
from Figures 2 and 3, and W B can be calculated by equation 



| = sin A 



(19) 



(17). 
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We obtain from equation (9) the circulatory velocity of a 
streamline point 

♦ (Y,X) = -.3R»l sin a th m J X s + Y 5 " 



= -2^1 sin a th j yg 

*1 



dZ 

W z | = 2 ^ V w .sin a th , (20) 



6ince 



7x 8 + V s = R (21) 



These equations contain no unknown quantities, since 
sin was already known. The resultant velocity is the sum 

of the parallel-flow velocity and the circulatory velocity. 
It should be noted, however, that the latter can have both pos- 
itive and negative values. This corresponds also to the purely 
logical conclusion, since the flow velocity on the lower side 
of the profile, due to opposite circulatory flow, is smaller 
than on the upper side where both flows have the 6ame direction 
and strengthen each other. 

The velocities W(Z), thus calculated for the flow around 
the cylinder must be still further corrected for the wing pro- 
file. This is possible, due to the following considerations. 
As already mentioned, the flow variations in the profile plane 
z, for the short intervals of time, are. proportional to the 
variations in the flow in the plane Z of the basic circle, when 
the variations in the Z values with respect to the z values are 
simultaneously considered. 
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dl = df dZ 
dz dZ dz 



or 



|w(z)| = 



W(Z) 
dz 
dZ 



w(z) 

1 - 1 



12 

(25) 

(26) 



because 



2=2— "|r- (functional equation of the 
Joukowsky function for con- 



formal transformation) 



(27) 



It then follows that 



dZ " x ' Z 2 '. 
The. divisor of equation (6) can be written 



1 




I 2 


1 - z 2 




z 



In the profile construction (Fig. 1) the line from the zero 
point to the outlying point on the streamline of the basic cir- 
cle equals Z, and the mirrored line to the point in the inner 

circle equals f/Z. Therefore, the line connecting the points 

7 2 

is the vectorial difference Z - ~. The quotient of the length 

a 

of the line Z divided by this difference is therefore the cor- 
rection factor by which the velocity on the circular profile 
must be multiplied, in order to obtain the correct velocity for 
the wing profile. 

The values of W and of the correction factor K for the 
wing profile are plotted separately in Figure 5, in order to 
show their characteristic^ courses in terms of their local rela- 
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tion to the circular profile. This representation shows very 
clearly the strong influence of the profile on the region near 
its edge and the general balancing of the velocity differences 
with increasing distance from the profile edge. In its funda- 
mental shape the curve of the correction factor forms the coun- 
terpart to the course of the velocity curves. The asymmetry of 
the profile shape is noticeable, however, this being the cause 
of the velocity variations as compared with the circular profile. 

TABLE III 



St re ami: 


me I 


Point No. 


1 


2 


3 


4 


5 


6 


7 


Wp/V^ 

Wz/Voo 


0.991 
0.119 


1.047 
0.132 


1.13 
0.143 


1.237 
0.154 


1.13 
0.143 


1.047 
0.132 


0.991 
0.119 


K 


1.110 
1.02 


1.179 
0.99 


1.273 
0.94 


1.391 
0-.865 


1.273 
0.89 


1.179 
0.94 


1.110 
0.98 




1.13 




1.2 


1.202 


1.13 


1.11 


1.09 








Streamline II 










0.88 
0.143 


0.954 
0.167 


1.144 

0.189 ; 


1.418 
0.206 


1.144 
0.189 


0.954 
0.167 


0.88 
0.143 


W/Voo 

K 


1.023 
1.1 


1.121 
1.05 


1.333 
0.97 


1.624 
0.79 


1.333 
0.87 


1.121 
0.97 


1.023 
1.06 




1.125 


1.179 


1.297 


1.284 


1.156 


1.109 


1.086 
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TABLE III (cont'd) 



|w| = IWl K 



Streamline III 



Point No. 


1 


2 


3 


4 


5 


6 


7 


Wp/V w 


0.75 
0.158 


0.592 
0.206 


0.912 , 
0.269 


1.826 
0.288 


0.608 
0.269 


0.592 
0.206 


0.75 
0.158 




0.908 


0. 798 

1. ob 


1.181 
1.16 


2.114 
0.67 


1.181 
u# by 


0.798 


0.908 

±9 <>0 




1.088 


1.089 


1.377 


1.413 


1.053 


l.ii 


1.14 








Streamline IV 








Wp/V^ 
Wz/V w 


0.942 
0.132 


1.002 
0.156 


1.15 
0.181 


1.396 
0.199 


1.15 
0.181 


1.002 
0.156 


0.942 
0.132 


K 


0.810 
1.05 


0.846 
1.01 


0.969 
0.91 


1.197 
0.75 


0.969 
0.90 


0.846 
1.03 


0.81 
1.08 




0.85 


0.857 


0.897 


0.898 


0. 876 


0.875 


0.872 








Streamline V 








w p /v w 


0.791 
0.154 


0.68 
0.202 


0.88 
0.27 


1.84 
0.291 


0.88 
0.27 


0.68 
0.202 


0.791 
0.154 


K 


0.637 
1.16 


•0.478 
1.29 


0.61 
1.26 


1.549 
0.56 


0.61 
1.39 


0.478 
1.52 


0.637 
1.24 




0.74 


0.619 


0.766 


0.875 


0.846 


0.726 


0.79 



Table III contains the velocities calculated for the five 
streamlines shown in Figure 5. The results are graphically rep- 
resented in Figure 6. The values liave already been introduced 
with the necessary distortion for the investigated wing profile 
at the corresponding locations, on the wing-profile streamlines, 
of the points coordinated with the circular streamlines. If 
cross sections are now made perpendicularly to the wing-profile 
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chord, it is then possible to plot the velocity drop for each 
section from the wing profile outward* The velocity on the con- 
tour of the wing profile itself can be easily calculated, since 
R = R x for this limiting case, so that equation (17) takes the 
form 



and equation (20) 



or 



= 2 V a sinA (22) 

= 2V a sin ath U 3 ) 
= 2 Vo, (sin A + sin ath) (24) 



For this limiting case Trefftz has developed a simplified method 
of calculation. (Reference 1,-page 18) ♦ 

As the abscissa of the last-mentioned diagram (Fig. 7), we 
chose the percentile distance of each streamline point from the 
profile edge, referred to the maximum profile thickness as unity 
and, as the ordinate, the velocity corresponding to each point. 
In this way we obtained, for any dimensions of the profile, a 
practical representation, which can also be used for determin- 
ing the velocity distribution for any propeller arrangement in 
the plane of any profile section. 

Analysis of the Results 

On comparing the increase in the velocities on the stream- 
lines I, II and III and the velocity drop of the streamlines 
from within outward, the regularity of the velocity distribution 
is obvious. The effect of the profile disturbance on the flow 
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increases from within outward, attaining its maximum value at the 
profile edge. The velocity increment in the direction. of a stream-* 
line is distributed in such a way that the velocity of a stream- 
line attains its highest value in the forward third of the pro- 
file and decreases from that point toward the trailing edge. 
The theoretically calculable maximum velocity is about 1.6 times 
the flow velocity* 

The effect of the profile on the flow is noticeable for some 
distance fore and aft of the profile. The velocity on the edge 
of the profile drops to zero at the points of full dynamic pres- 
sure. Its effect can still be observed in the velocity of the 
streamlines near the profile edge, since there is a considerable 
velocity drop in these streamlines fore and aft of the profile. 
Consequently there is, in the streamlines fore and. aft of this 
point a further velocity increase, before this velocity drops, 
at the corresponding distance, to its normal value of %>. In 
more remote streamlines a certain, equalization has already taken 
place, so that a constant velocity drop is observable with in- 
creasing distance from the profile. Like phenomena axe observ- 
able in streamlines IV and V on the lower ,side v of the profile. 

If the individual planes of the profile sections are exam- 
ined in the light of these results, there is to be expected, for 
propellers mounted above the profile, a relatively uniform ve- 
locity distribution in the intersecting plane a (Fig. 7) in 
front of the profile. The plane b at the height of the lead- 
ing edge of the profile is also favorable for propellers, provid- 
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ed the blade tips do not come too close to the profile edge* 
In plane c, on the contrary, great velocity differences and 
corresponding unequal deflections of the propeller-disk surfaces 
are to be expected. Conditions similar to those for c obtain 
for planes d to g, but the conditions gradually improve to- 
ward the trailing edge of the wing. Near the profile edge the 
plane h shows a certain velocity drop, which probably has no 
effect, however, at the customary minimum distance. The trailing 
edge of the profile offers very favorable velocity relations. 
The plane k shows, at about the height of. the trailing edge, 
a small velocity increase for this portion of the field, the 
efficiency of which depends on the height of the propeller. 

The conditions are similar on both sides of the profile. 
The only exceptions are the cross sections a and i, whose ve- 
locity fields are subject to certain irregularities* The planes 
c, g, ojid h likewise exhibit certain velocity differences, but 
they are not nearly so great as on the upper side. 

The calculated order of magnitude of the velocity variations 
must naturally be adopted with caution and should first be veri- 
fied by practical experiments, since practical conditions often 
differ from the ideal theoretical conditions and produce effects 
which cannot be taken care of mathematically. It might be con- 
cluded from the results of the Gottingen pressure measurements 
on wing profiles that the tendency of the curve is correct, since 
the distribution of the pressure along the profile chord was 
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there found to be similar to that calculated here for the veloc- 
ity distribution, and since the pressure distribution, as kinetic 
energy, is a function of the velocity. 

A knowledge of the profile effect on the flow velocity is 
also of some importance for the application of dynamic-pressure 
meters, because it enables one to estimate the minimum distance 
from the wing for measuring the least possible affected flow ve- 
locity. 
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